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Abstract 

We present results of a search for the rare flavor-changing neutral-current decay B — ► 7r^ + £~, 
based on a data sample corresponding to 209 fb _1 of integrated luminosity collected with the BABAR 
detector at the PEP-II B Factory. We reconstruct the four exclusive B decay modes B + — * ir + £ + £~ 
and B° — ► Tt°£ + £~ , where £ is either an e or \i. We find no evidence for a signal, and we obtain the 
upper limit at 90% confidence level on the lepton-flavor-averaged branching fraction to be 

B{B + -» n + £ + £~) = 2 x ^B(B° -» vr°£ + r) < 7.9 x 10~ 8 . 

We also obtain an upper limit at 90% confidence level on the lepton-flavor-violating decay B — ► 7re/x 
of 

S(-B vreu) < 9.8 x 10~ 8 . 

Submitted to the 33 rd International Conference on High-Energy Physics, ICHEP 06, 
26 July — 2 August 2006, Moscow, Russia. 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 
Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, R. Barate, M. Bona, D. Boutigny, F. Couderc, Y. Karyotakis, J. P. Lccs, V. Poireau, 

V. Tisserand, A. Zghichc 

Laboratoire de Physique des Particules, IN2P3/CNRS et Universite de Savoie, F-74941 Annecy-Le-Vieux, 

France 

E. Grauges 

Universitat de Barcelona, Facultat de Fisica, Departament ECM, E- 08028 Barcelona, Spain 

A. Palano 

Universita di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigcn, I. Ofte, B. Stugu 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

G. S. Abrams, M. Battaglia, D. N. Brown, J. Button-Shafcr, R. N. Cahn, E. Charles, M. S. Gill, 
Y. Groysman, R. G. Jacobsen, J. A. Kadyk, L. T. Kerth, Yu. G. Kolomensky, G. Kukartsev, G. Lynch, 
L. M. Mir, T. J. Orimoto, M. Pripstcin, N. A. Roe, M. T. Ronan, W. A. Wenzel 

Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94-720, USA 

P. del Amo Sanchez, M. Barrett, K. E. Ford, A. J. Hart, T. J. Harrison, C. M. Hawkes, S. E. Morgan, 

A. T. Watson 

University of Birmingham, Birmingham, B15 2TT, United Kingdom 

T. Held, H. Koch, B. Lcwandowski, M. Pclizaeus, K. Peters, T. Schroeder, M. Stcinke 
Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 

J. T. Boyd, J. P. Burke, W. N. Cottingham, D. Walker 
University of Bristol, Bristol BS8 1TL, United Kingdom 

D. J. Asgeirsson, T. Cuhadar-Donszelmann, B. G. Fulsom, C. Hearty, N. S. Knecht, T. S. Mattison, 

J. A. McKenna 

University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 

A. Khan, P. Kyberd, M. Saleem, D. J. Sherwood, L. Teodorescu 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

V. E. Blinov, A. D. Bukin, V. P. Druzhinin, V. B. Golubev, A. P. Onuchin, S. I. Serednyakov, 
Yu. I. Skovpcn, E. P. Solodov, K. Yu Todyshev 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. S. Best, M. Bondioli, M. Bruinsma, M. Chao, S. Curry, I. Eschrich, D. Kirkby, A. J. Lankford, P. Lund, 
M. Mandelkcrn, R. K. Mommsen, W. Roethcl, D. P. Stoker 

University of California at Irvine, Irvine, California 92697, USA 

S. Abachi, C. Buchanan 
University of California at Los Angeles, Los Angeles, California 90024, USA 



2 



S. D. Foulkes, J. W. Gary, O. Long, B. C. Shen, K. Wang, L. Zhang 
University of California at Riverside, Riverside, California 92521, USA 

H. K. Hadavand, E. J. Hill, H. P. Paar, S. Rahatlou, V. Sharma 
University of California at San Diego, La Jolla, California 92093, USA 

J. W. Berryhill, C. Campagnari, A. Cunha, B. Dahmes, T. M. Hong, D. Kovalskyi, J. D. Richman 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

T. W. Beck, A. M. Eisner, C. J. Flacco, C. A. Heusch, J. Kroseberg, W. S. Lockman, G. Nesom, T. Schalk, 
B. A. Schumm, A. Seiden, P. Spradlin, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

J. Albert, E. Chen, A. Dvorctskii, F. Fang, D. G. Hitlin, I. Narsky, T. Piatenko, F. C. Porter, A. Ryd, 

A. Samuel 

California Institute of Technology, Pasadena, California 91125, USA 

G. Mancinelli, B. T. Meadows, K. Mishra, M. D. Sokoloff 

University of Cincinnati, Cincinnati, Ohio 45221, USA 

F. Blanc, P. C. Bloom, S. Chen, W. T. Ford, J. F. Hirschauer, A. Kreisel, M. Nagel, U. Nauenberg, 
A. Olivas, W. O. Ruddick, J. G. Smith, K. A. Ulmer, S. R. Wagner, J. Zhang 

University of Colorado, Boulder, Colorado 80309, USA 

A. Chen, E. A. Eckhart, A. Soffcr, W. H. Toki, R. J. Wilson, F. Winklmeier, Q. Zcng 
Colorado State University, Fort Collins, Colorado 80523, USA 

D. D. Altenburg, E. Feltresi, A. Hauke, H. Jasper, J. Merkel, A. Petzold, B. Spaan 
Universitat Dortmund, Institut fur Physik, D-44221 Dortmund, Germany 

T. Brandt, V. Klose, H. M. Lacker, W. F. Mader, R. Nogowski, J. Schubert, K. R. Schubert, R. Schwierz, 

J. E. Sundcrmann, A. Volk 

Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonneaud, E. Latour, Ch. Thicbaux, M. Verderi 
Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Poly technique, F-91128 Palaiseau, France 

P. J. Clark, W. Gradl, F. Muhcim, S. Playfer, A. I. Robertson, Y. Xic 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinctto, E. Luppi, M. Negrini, A. Pctrclla, 

L. Piemontese, E. Prencipe 

Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, S. Pacetti, P. Patteri, 

I. M. Pcruzzi, 1 M. Piccolo, M. Rama, A. Zallo 

Laboratori Nazionali di Frascati dell' INFN, 1-00044 Frascati, Italy 

1 Also with Universita di Perugia, Dipartimento di Fisica, Perugia, Italy 



3 



A. Buzzo, R. Capra, R. Contri, M. Lo Vetere, M. M. Macri, M. R. Monge, S. Passaggio, C. Patrignani, 

E. Robutti, A. Santroni, S. Tosi 

Universitd di Genova, Dipartimento di Fisica and INFN, 1-16146 Genova, Italy 

G. Brandenburg, K. S. Chaisanguanthum, M. Morii, J. Wu 
Harvard University, Cambridge, Massachusetts 02138, USA 

R. S. Dubitzky, J. Marks, S. Schenk, U. Uwer 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

D. J. Bard, W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, R. L. Flack, J. A. Nash, 

M. B. Nikolich, W. Panduro Vazquez 

Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Bchera, X. Chai, M. J. Charles, U. Mallik, N. T. Meyer, V. Ziegler 
University of Iowa, Iowa City, Iowa 52242, USA 

J. Cochran, H. B. Crawley, L. Dong, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin 
Iowa State University, Ames, Iowa 50011-3160, USA 

A. V. Gritsan 

Johns Hopkins University, Baltimore, Maryland 21218, USA 

A. G. Denig, M. Fritsch, G. Schott 
Universitdt Karlsruhe, Institut fur Experimentelle Kernphysik, D- 76021 Karlsruhe, Germany 

N. Arnaud, M. Davier, G. Grosdidier, A. Hocker, F. Le Diberder, V. Lepeltier, A. M. Lutz, A. Oyanguren, 
S. Pruvot, S. Rodier, P. Roudeau, M. H. Schune, A. Stocchi, W. F. Wang, G. Wormser 

Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, Centre Scientifique 

d'Orsay, B.P. 34, F-91898 ORSAY Cedex, France 

C. H. Cheng, D. J. Langc, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

C. A. Chavez, I. J. Forster, J. R. Fry, E. Gabathuler, R. Garnet, K. A. George, D. E. Hutchcroft, 

D. J. Payne, K. C. Schoficld, C. Touramanis 

University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bcvan, F. Di Lodovico, W. Menges, R. Sacco 
Queen Mary, University of London, El 4NS, United Kingdom 

G. Cowan, H. U. Flaecher, D. A. Hopkins, P. S. Jackson, T. R. McMahon, S. Ricciardi, F. Salvatore, 

A. C. Wren 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United 

Kingdom 

D. N. Brown, C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 



4 



J. Allison, N. R. Barlow, R. J. Barlow, Y. M. Chia, C. L. Edgar, G. D. Laffcrty, M. T. Naisbit, 

J. C. Williams, J. I. Yi 

University of Manchester, Manchester MIS 9PL, United Kingdom 

C. Chen, W. D. Hulsbcrgen, A. Jawahery, C. K. Lae, D. A. Roberts, G. Simi 
University of Maryland, College Park, Maryland 20742, USA 

G. Blaylock, C. Dallapiccola, S. S. Hertzbach, X. Li, T. B. Moore, S. Saremi, H. Staenglc 
University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, G. Sciolla, S. J. Sekula, M. Spitznagel, F. Taylor, R. K. Yamamoto 

Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, 

USA 

H. Kim, S. E. Mclachlin, P. M. Patel, S. H. Robertson 
McGill University, Montreal, Quebec, Canada H3A 2T8 

A. Lazzaro, V. Lombardo, F. Palombo 
Universita di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Krocgcr, D. A. Sanders, D. J. Summers, 

H. W. Zhao 

University of Mississippi, University, Mississippi 38611, USA 

S. Brunet, D. Cote, M. Simard, P. Taras, F. B. Viaud 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3.11 

H. Nicholson 

Mount Holyoke College, South Hadley, Massachusetts 01015, USA 

N. Cavallo, 2 G. De Nardo, F. Fabozzi, 3 C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, 

C. Sciacca 

Universita di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. A. Baak, G. Raven, H. L. Snoek 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The 

Netherlands 

C. P. Jessop, J. M. LoSecco 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

T. Allmcndingcr, G. Bcnclli, L. A. Corwin, K. K. Gan, K. Honschcid, D. Hufnagcl, P. D. Jackson, 
H. Kagan, R. Kass, A. M. Rahimi, J. J. Regensburger, R. Ter-Antonyan, Q. K. Wong 

Ohio State University, Columbus, Ohio 43210, USA 

N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, M. Lu, R. Rahmat, N. B. Sinev, D. Strom, 

J. Strubc, E. Torrence 

University of Oregon, Eugene, Oregon 91403, USA 

2 Also with Universita della Basilicata, Potenza, Italy 
3 Also with Universita della Basilicata, Potenza, Italy 



5 



A. Gaz, M. Margoni, M. Morandin, A. Pompili, M. Posocco, M. Rotondo, F. Simonetto, R. Stroili, C. Voci 
Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, L. Del Buono, Ch. de la Vaissiere, O. Hamon, 
B. L. Hartfiel, M. J. J. John, Ph. Leruste, J. Malcles, J. Ocariz, L. Roos, G. Therin 

Laboratoire de Physique Nucleaire et de Hautes Energies, IN2P3/CNRS, Universite Pierre et Marie 
Curie-Paris6, Universite Denis Diderot- Paris!, F-75252 Paris, France 

L. Gladney, J. Panetta 
University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

M. Biasini, R. Covarelli 
Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

C. Angelini, G. Batignani, S. Bettarini, F. Bucci, G. Caldcrini, M. Carpinclli, R. Ccnci, F. Forti, 
M. A. Giorgi, A. Lusiani, G. Marchiori, M. A. Mazur, M. Morganti, N. Ncri, E. Paoloni, G. Rizzo, 

J. J. Walsh 

Universita di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

M. Haire, D. Judd, D. E. Wagoner 
Prairie View A&M University, Prairie View, Texas 77446, USA 

J. Biesiada, N. Daniclson, P. Elmer, Y. P. Lau, C. Lu, J. Olsen, A. J. S. Smith, A. V. Telnov 
Princeton University, Princeton, New Jersey 08544, USA 

F. Bellini, G. Cavoto, A. D'Orazio, D. del Re, E. Di Marco, R. Faccini, F. Ferrarotto, F. Ferroni, 
M. Gaspero, L. Li Gioi, M. A. Mazzoni, S. Morganti, G. Piredda, F. Polci, F. Safai Tehrani, C. Voena 

Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

M. Ebert, H. Schroder, R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, N. De Groot, B. Franck, E. O. Olaiya, F. F. Wilson 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, G. Hamel de Monchenault, W. Kozanecki, M. Legendre, 

G. Vasseur, Ch. Yeche, M. Zito 

DSM/Dapnia, CEA/Saclay, F- 91 191 Gif-sur-Yvette, France 

X. R. Chen, H. Liu, W. Park, M. V. Purohit, J. R. Wilson 
University of South Carolina, Columbia, South Carolina 29208, USA 

M. T. Allen, D. Aston, R. Bartoldus, P. Bcchtle, N. Berger, R. Claus, J. P. Coleman, M. R. Convery, 

M. Cristinziani, J. C. Dingfelder, J. Dorfan, G. P. Dubois-Felsmann, D. Dujmic, W. Dunwoodie, 
R. C. Field, T. Glanzman, S. J. Gowdy, M. T. Graham, P. Grenier, 4 V. Halyo, C. Hast, T. Hryn'ova, 

W. R. Innes, M. H. Kelsey, P. Kim, D. W. G. S. Leith, S. Li, S. Luitz, V. Luth, H. L. Lynch, 
D. B. MacFarlane, H. Marsiske, R. Messner, D. R. Muller, C. P. O'Grady, V. E. Ozcan, A. Perazzo, 
M. Perl, T. Pulliam, B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwiening, 
A. Snyder, J. Stelzer, D. Su, M. K. Sullivan, K. Suzuki, S. K. Swain, J. M. Thompson, J. Va'vra, N. van 

4 Also at Laboratoire de Physique Corpusculaire, Clermont-Ferrand, France 



6 



Bakel, M. Weaver, A. J. R. Weinstein, W. J. Wisniewski, M. Wittgen, D. H. Wright, A. K. Yarritu, K. Yi, 

C. C. Young 

Stanford Linear Accelerator Center, Stanford, California 94-309, USA 

P. R. Burchat, A. J. Edwards, S. A. Majewski, B. A. Petersen, C. Roat, L. Wilden 
Stanford University, Stanford, California 94305-4060, USA 

S. Ahmed, M. S. Alam, R. Bula, J. A. Ernst, V. Jain, B. Pan, M. A. Saeed, F. R. Wappler, S. B. Zain 
State University of New York, Albany, New York 12222, USA 

W. Bugg, M. Krishnamurthy, S. M. Spanicr 

University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. Satpathy, C. J. Schilling, R. F. Schwitters 
University of Texas at Austin, Austin, Texas 18112, USA 

J. M. Izen, X. C. Lou, S. Ye 

University of Texas at Dallas, Richardson, Texas 15083, USA 

F. Bianchi, F. Gallo, D. Gamba 
Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

M. Bomben, L. Bosisio, C. Cartaro, F. Cossutti, G. Delia Ricca, S. Dittongo, L. Lanccri, L. Vitalc 
Universitd di Trieste, Dipartimento di Fisica and INFN, 1-34121 Trieste, Italy 

V. Azzolini, N. Lopez-March, F. Martinez- Vidal 
IFIC, Universitat de Valencia- CSIC, E-46011 Valencia, Spain 

Sw. Banerjee, B. Bhuyan, C. M. Brown, D. Fortin, K. Hamano, R. Kowalewski, I. M. Nugent, J. M. Roney, 

R. J. Sobie 

University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

J. J. Back, P. F. Harrison, T. E. Latham, G. B. Mohanty, M. Pappagallo 
Department of Physics, University of Warwick, Coventry CV4 1AL, United Kingdom 

H. R. Band, X. Chen, B. Cheng, S. Dasu, M. Datta, K. T. Flood, J. J. Hollar, P. E. Kuttcr, B. Mcllado, 
A. Mihalyi, Y. Pan, M. Picrini, R. Prepost, S. L. Wu, Z. Yu 

University of Wisconsin, Madison, Wisconsin 53106, USA 

H. Neal 

Yale University, New Haven, Connecticut 06511, USA 



7 



1 Introduction 



The decays B — > 7r£ + £ _ , where is either an e + e~ or p + p~ pair, are the simplest manifestations 
of b — » d£ + £~ flavor-changing neutral currents (FCNC). In the Standard Model (SM), these decays 
are forbidden at tree level and can only occur at greatly suppressed rates through higher-order 
processes. At lowest order in the electroweak couplings, three amplitudes contribute: (1) a photon 
penguin, (2) a Z penguin, and (3) a box diagram (Fig. Since these decays proceed 

via weakly-interacting particles with virtual energies near the electroweak scale, they provide a 
promising means to search for effects from new flavor-changing interactions. Such effects are pre- 
dicted in a wide variety of models, usually in the context of the similar FCNC b — > s£ + £~ [1-5]. 
If there exists non-trivial flavor violation in the new interactions, b — * d£ + £~ can also exhibit large 
observable effects, independent of the experimental constraints on b — > s£ + £~ [6,7]. 

In the Standard Model, the b — > d FCNC decay rate is suppressed relative to b — > s by 
a ratio of Cabibbo-Kobayashi-Maskawa matrix elements (V^/V^l 2 , making them especially rare 
processes which challenge the sensitivity of the B factories. The chief experimental constraint on 
AB = 1 b — » d FCNC decays comes from recent measurements of the exclusive decays B —* pj and 
B — > ujj [8], which constrain the photon penguin decay rate to agree with SM predictions with an 
experimental precision of 30%. The additional b — > d electroweak amplitudes present in b — > d£ + £~ 
are weakly constrained experimentally from older measurements [9]. A recent SM prediction [10] 
for the B + — * n + £ + £~ branching fraction is based on an operator product expansion combined with 
B — » 7r form factor predictions from light-cone QCD sum rules; the predicted branching fraction 
is 3.3 x 10 -8 with an uncertainty of about 30%. The B° — ► ir°£ + £~ branching fraction is related 
to the charged mode by isospin symmetry; accounting for the difference in charged and neutral B 
meson lifetimes, the branching fraction is given by \ x t b o/t b + x B(B+ -> 7T+ £+£-). 




q q q q 

Figure 1: Standard Model diagrams for the decays B — > ir£ + £~. For the photon or Z penguin 
diagrams on the left, boson emission can occur on any of the b, t, c, u, d, or W lines. 

We also perform a search for the lepton-flavor violating decay B — > ire^fi^ which can occur 
in some models beyond the Standard Model, such as theories involving leptoquarks [5]. Earlier 
searches [11] for these modes have yielded no such events. 
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2 The BABAR Detector and Dataset 



The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring 
located at the Stanford Linear Accelerator Center. The data sample comprises 208.9 fb _1 recorded 
near the peak of the T(4S) resonance, yielding (230.1 ±2.5) x 10 6 BB decays, and an off-resonance 
sample of 21.5 fb -1 which is used to study continuum background. 

The BABAR detector is described in detail elsewhere [12]. Its most important capabilities for this 
study are charged-particle tracking and momentum measurement, charged tt/K separation, and lep- 
ton identification. Charged particle tracking is provided by a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH). The DIRC, a Cherenkov ring- imaging particle-identification 
system, is used to identify charged kaons and pions. Electrons are identified using the electromag- 
netic calorimeter (EMC), which comprises 6580 thallium-doped Csl crystals. These systems are 
mounted inside a 1.5 T solenoidal superconducting magnet. Muons are identified in the instru- 
mented flux return (IFR), in which resistive plate chambers are interleaved with the iron plates of 
the magnet flux return. 

3 Analysis Method 
3.1 Event Selection 

We select events that include two oppositely charged leptons (e + e~, fi + p~ and e ± ^x ¥ ) as well as 
a pion (either a 7r =l= track or a ir° meson decaying to 77). The analysis faces the experimental 
challenge of isolating a very rare B-meson decay from background due to random combinations of 
particles in generic e + e~ — > T(4S) — ► BB and continuum e + e~ — ► qq events (qq being a pair of 
u,d,s, or c quarks), and there is a great abundance of pions and leptons in all of these events. We use 
stringent particle identification and a multivariate algorithm to reduce this background. Another 
type of background comes from events that have signal-like features. The main contributions to 
such backgrounds are from B-mesons decaying to charmonium states (e.g. B — ► J/tp it where 
the J/ip decays to two leptons) or B-mesons decaying to three hadrons (e.g. B — > Dtt with 
D — > Kir) where two of the hadrons pass the muon selection. In addition, there is background 
from the b — ► s£ + £~ penguin decay B — > K £ + £~ , which has a branching fraction about one order 
of magnitude larger than expected for B — > ir £ + £~ [6, 10]. These events have a reconstructed B 
energy lower than expected for signal events and are also very efficiently reduced by information 
from the DIRC subdetector. 

Electrons are required to have a momentum p > 0.3 GeV/c and are identified combining infor- 
mation from the EMC, DIRC, and DCH. Photons with E > 30 MeV that lie within an angular 
region of 50 mrad in the azimuth angle <p and 35 mrad in the polar angle 9 of the initial electron track 
direction are combined with electron candidates in order to recover energy lost by bremsstrahlung. 
We suppress backgrounds due to photon conversions in the B — > tt e + e~ channels by removing 
e + e~ pairs with invariant mass less than 30 MeV/c 2 . Muons with momentum p > 0.7 GeV/c are 
identified with a neural network algorithm using information from the IFR and the EMC. 

The lepton pairs are combined with a charged or neutral pion candidate to reconstruct a B 
meson. Charged pion candidates are tracks with specific ionization (dE/dx) and Cherenkov angle 
consistent with a pion. A likelihood algorithm based on this information reduces fake rates from 
kaons down to about 1 — 5% depending on the track momentum. Neutral pion candidates are 
identified as pairs of neutral-energy deposits in the EMC, each with an energy greater than 50 
MeV in the laboratory frame. The invariant mass of the pair is calculated under the assumption 
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that the photons originate from the £ + £ vertex, and is required to satisfy 115 < m 77 < 150 
MeV/c 2 . 

Correctly reconstructed B signal decays produce narrow peaks in the distributions of two kine- 
matic variables, which can be used to extract the signal and background yields. For a candidate 
system of B daughter particles with masses rrii and three-momenta p* in the T(4S') center-of-mass 

(CM) frame, we define m ES = \j E^ 2 - |EiP*| 2 and AE = J2i y^i + Pf ~ K> where E b is the 
beam energy in the CM frame. For signal events, the tties distribution peaks at the B meson mass 
and has a width <7 mBS ~ 2.5 MeV/c 2 , and the AE distribution peaks near zero, with a typical width 
cab ~ 20 MeV. The mean and width of the tues and AE distributions are determined separately 
for e + e~ and modes using data control samples. For events reconstructed as e^ix^ we assume 

the same width and mean as for the e + e~ modes. 

We define two kinematic regions in terms of th-es and AE for signal extraction purposes. The 
signal region is defined to be within 2a of the expected mean of the peak in mES and AE. The 
values are given in Tabled A 2a signal region width was found to be close to optimal for both the 
Sfy/S + B and S/y/B figure of merit, where S and B are the number of signal and background 
events expected based on Monte Carlo simulations. 

Table 1: Boundary values defining the signal region for each B — > ir£ + £~ mode. The boundaries 
used in the B + — > 7r + e / u and B° — ► 7r°e/i modes are the same as for B + — » 7r + e + e~ and B° — > 
7r°e + e~ modes, respectively. 



mode 


mES low 
[GeV/c 2 ] 


m E s high 
[GeV/c 2 ] 


AE low 

[MeV] 


AE high 
[MeV] 


B + -> 7r + e + e" 


5.2748 


5.2847 


-53.6 


37.4 


B° -> 7r°e+e- 


5.2767 


5.2839 


-115.0 


82.5 


B + -» 7T + /i+^- 


5.2749 


5.2847 


-42.0 


35.0 


B° -> ttV + M" 


5.2764 


5.2836 


-87.4 


68.0 



We blind ourselves from inspecting the signal region, and also a broader region in the mES"^^ 
plane, until all selection criteria are optimized and control samples are checked. The broader region 
is used to determine the number of background events from an unbinned maximum likelihood fit 
directly from the data. This fit region is defined as 5.2GeV/c 2 < m-Es and \ AE\ < 0.25 GeV. 

Backgrounds arise from three main sources: random combinations of particles from qq and 
t + t~ events produced in the continuum, random combinations of particles from T(4S) — > BB 
decays, and B decays to topologies similar to the signal modes. The first two, "combinatorial", 
backgrounds typically arise from pairs of semileptonic decays and produce broad distributions in 
mES and AE compared to the signal. The third source arises from modes which have shapes similar 
to the signal, such as B — > J/ifm, with J ftp — > £ + £~ , or B — > Ktttt, with kaons or pions misidentified 
as muons. All selection criteria are optimized with GEANT4 [13] simulated data, independent data 
control samples or with data samples outside the full fit region. 

We suppress combinatorial background from continuum processes by selecting events with high 
values in a Fisher discriminant [14] constructed as a linear combination of variables with coefficients 
optimized to distinguish between signal and background. The variables (defined in the CM frame) 
are (1) the ratio of second- to zeroth-order Fox- Wolfram moments [15] for the event, computed 
using all charged tracks and neutral energy clusters; (2) the absolute value of the angle between the 
thrust axis of the B candidate and that of the remaining particles in the event; (3) the production 
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angle 9b of the B candidate with respect to the beam axis; and (4) the ratio of second- to zeroth- 
order Legendre moments [16] for the event, computed using all charged tracks and neutral energy 
clusters. These variables exploit the difference between the relatively spherical track distribution 
in BB events and the jet-like structure of continuum events. 

We suppress combinatorial backgrounds from BB events by selecting events with high values 
in a likelihood function constructed from (1) the missing energy of the event, computed from all 
charged tracks and neutral energy clusters; (2) the vertex fit probability of all tracks from the 
B candidate; (3) the vertex fit probability of the two leptons; and (4) the angle 9b of the B 
candidate with respect to the beam axis in the direction of the negative electron beam. Missing 
energy provides the strongest suppression of combinatorial BB background events, which typically 
contain energetic neutrinos from at least two semileptonic B or D meson decays. 

Both the continuum and BB suppression discriminants are constructed separately for each e + e~ 
and fi~ mode. Parameters of the Fisher discriminant have been determined from simulated sig- 
nal events and a data control sample selected below the T(4S) energy. The likelihood function 
parameters have been determined from simulated signal events and simulated generic BB events 
to separate signal events from other types of BB backgrounds. For each reconstruction mode, the 
required minimum values of the Fisher and likelihood functions are optimized simultaneously by 
maximizing S/yfS + B, where S is the expected signal yield in the signal region, and B is the 
expected background yield in the signal region based on a combination of generic BB and contin- 
uum Monte Carlo simulated events. The background yields are extracted from a two-dimensional, 
unbinned, maximum-likelihood fit to the simulated events in the fit region. A branching fraction 
of B(B — > tt££) of 3.3 x 10~ 8 [10] is assumed as the signal estimate. The optimal selection criteria 
do not change significantly if we alternatively optimize for S/yB. For the ire^^ sample we use 
the same parameters and optimal selection points as found for the e + e~ modes. 

The most prominent backgrounds that peak in the signal region are B decays to charmonium: 
B — > J/i/j 7r (with J/i[) — » The more abundant B — > J '/ip K decay also contributes as 

background to B —* n l + t~ . This latter background is strongly suppressed by particle identi- 
fication requirements, and such events typically have a reconstructed B energy which is lower 
than expected for our signal mode, so that most of these events are outside the signal region for 
B — > it t + £~ . Similarly, B decays to i{)(2S)ir(K) final states contribute with a peaking back- 
ground component. We exclude dilepton pairs with dilepton invariant mass mg+£- consistent with 
the J/V> mass (2.90 < m e+e - < 3.20 GeV/c 2 and 3.00 < m^ +fl - < 3.20 GeV/c 2 ) or with the 
ip(2S) mass (3.60 < rri£+£- < 3.75 GeV/c 2 ). This veto is also applied to m e + e - computed without 
bremsstrahlung photon recovery. When a lepton radiates or is mismeasured, mg+£- can shift away 
from the charmonium mass, while AE shifts in a correlated manner. Therefore, we apply an ad- 
ditional veto with rri£+£- veto regions defined as a linear function of AE. By extending the veto 
region this way, we remove nearly all charmonium events from the fit region and simplify the shape 
of the background PDF in the fit region. For e^/x^ the charmonium vetoes are defined to be the 
same as for the e + e~ modes. 

The charmonium events removed by these vetoes are otherwise kinematically similar to signal 
events. These events therefore serve as copious control samples for studying signal shapes, selec- 
tion efficiencies, and systematic errors. Figure |2 shows the distributions of the two background 
rejection variables for the B + —* J/ip K + control sample. This sample is selected with identical 
requirements as our signal sample, except the J ftp veto has been reversed and the hadron track 
passes kaon identification. The figure also shows the distributions for simulated B + — ► J/iJj; K + 
events and the simulations agree well with the data. We use this sample to measure the efficiency 



11 



and bound systematic uncertainties of lepton identification, Fisher discriminant, and B likelihood 
requirements. 

In addition to studying the efficiencies using the charmonium control sample, we measure the 
branching fractions of B(B — > J/ip it) using B — ► tt£ + £~ candidates which fail the charmonium 
veto. The measured branching fractions are consistent with the world averages for these modes [17]. 
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Figure 2: Distribution of the Fisher and likelihood variables in the J/tpK + control sample. The 
black points are BABAR data events and the histograms show comparison with simulated events, 
with statistical and systematic uncertainties shown in the dark and light gray bands, respectively. 
The vertical lines indicate the optimal cut values for each variable. Our event selection corresponds 
to events to the right of the vertical line. 



In modes with two muons in the final state, where the probability for a hadron to be misidentified 
as a muon can be as high as a few percent, background from hadronic -B-decays like B — ► Dir 
with D — ► 7T7T or D — > Kir is significant. These events are suppressed by vetoing events where 
the B — > it £ + l~ kinematics are consistent with those of a hadronic D decay. We evaluate the 
mg+g- and m n £ invariant masses with tt or K mass hypotheses on the tracks in all combinations 
corresponding to known D-meson decay modes; we veto those events which have two oppositely- 
charged tracks with an invariant mass within the range 1.84 — 1.89 GeV/c 2 , or events which have 
a neutral pion and charged track combination with an invariant mass within the range 1.79 — 1.94 
GeV/c 2 . 

If more than one charged B-meson candidate remain after all selection requirements have been 
applied, we select the candidate for which the pion has the greatest number of hits in the SVT 
tracking detector. If more than one neutral 5-meson candidate remain, we select the first one 
that appears in our data sample (i.e. a random candidate); the efficiency of this choice is nearly 
identical to choosing based on the reconstructed 7r° mass. These criteria are chosen in favor of 
methods using kinematic information to avoid any bias due to possible correlations with the signal 
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extraction variables toes an d AE. More than one B candidate per event occurs in about 13- 
15% of charged--B signal events and 35-60% of neutral-S signal events, and typically the multiple 
candidates are due to more than one possible 7r candidate. For the e fi^ modes more than one B 
candidate occurs in about 5% of events. 

3.2 Signal Extraction and Background Estimate 

The signal is extracted by counting events in the signal region. An upper limit on the branching 
fraction B(B — ► tt £ + £~) is determined from the observed number of events and the expected 
number of background events in this region. We determine the background expectation before 
inspecting the signal region. 

The number of combinatorial background events is extracted by performing a two-dimensional 
unbinned maximum- likelihood fit to toes and AE in the region toes < 5.2724 GeV/c 2 and | A.£7| < 
0.25 GeV. The background probability distribution function (PDF) is modeled as the product of an 
ARGUS function [18] for toes an d an exponential function for AE. The slopes and normalization 
of these functions are floating in the fit, while the endpoint of the ARGUS function is fixed to 
5.290 GeV/c 2 . The PDF is extrapolated into the signal region (up to toes < 5.290 GeV/c 2 ) and 
the expected number of background events in the signal region is obtained from integrating the 
PDF over the signal region. The expected number of combinatorial background events is given in 
Table 

We estimate the remaining hadronic B background from a data control sample of mainly B — > 
tttttt and B — > Ktttt events. In the signal sample these events are highly suppressed by the lepton 
identification criteria. For the hadronic control sample we select events where the charged tracks fail 
electron and muon identification, except that for data reduction purposes we require that one track 
(and only one) pass a looser muon identification criteria which still has a high hadron fake rate. 
This selects a sample of predominantly hadronic events. These hadronic events are further weighted 
by their probability to pass the signal-selection muon identification, by weighting each track with 
a probability depending on its particle type hypothesis, momentum and direction. These track-by- 
track probabilities are determined from data control samples specifically selected to study particle 
identification efficiencies. Looking at events with a AE consistent with a correctly reconstructed 
5-meson, we fit the weighted toes distribution for these events to obtain the expected number of 
hadronic peaking background events in the signal region. We estimate a background of 0.06 ± 0.05 
events in the ir + fi + fi" channel and 0.11 ± 0.04 events in the tt h + h~ channel, where the dominant 
uncertainty is the statistical uncertainty of the fit. We expect no such peaking background in the 
electron modes or in the e\x modes. 

£>-meson decays to K£ + £~ and p£ + £~ in the final state are the only other major peaking 
background components expected. The former shift toward lower reconstructed B energy and are 
mostly outside the 2a signal region. This is the sole largest peaking background component for the 
7r + e + e~ mode. The latter contribute even less, since the reconstructed B in these cases is missing 
a pion. These estimates are based on high-statistics samples of simulated events. In total, we 
estimate the number of peaking background events from leptonic events to range from 0.07 ± 0.02 
events for the 7r + e + e _ mode down to 0.02 ±0.01 events for the 7r°/i + ^~ mode. No leptonic peaking 
background is expected in the e\i modes. 

After the vetoes of B — > J/tpir(K) and B —* ip(2S) ir(K) decays, no remaining background 
from these are expected, as estimated from Monte Carlo simulations. 

Table [2] summarizes the number of background events expected in the signal region and the 
systematic uncertainties associated with these estimates. 
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Table 2: Estimated number of background events and their associated systematic uncertainties. 
The mss-AE fit uncertainties are evaluated by varying the parameters of the fit by ±lcx, and the 
mgs-Ai? correlations and AE shape uncertainties are evaluated from using alternative probability 
density functions. Estimates and systematics from the peaking background are based on Monte 
Carlo and control sample studies. 



Systematic 


7r + e + e 


7T U e+e- 


7T + fl + fl 


ttW 


7r + efj, 


7r u e/i 


m-Es-AE fit 


0.89 ± 0.31 


0.43 ±0.21 


0.86 ± 0.24 


0.22 ±0.18 


1.48 ±0.33 


1.13 ±0.40 


tues-AE correlations 


±0.02 


±0.03 


±0.06 


±0.03 


±0.17 


±0.05 


AE shape 


±0.03 


±0.03 


±0.15 


±0.02 


±0.31 


±0.24 


Leptonic peaking bkg. 


0.07 ±0.02 


0.03 ±0.01 


0.04 ±0.01 


0.02 ± 0.01 


0.00 ±0.00 


0.00 ±0.00 


Hadronic peaking bkg. 


0.00 ± 0.00 


0.00 ±0.00 


0.06 ±0.05 


0.11 ± 0.04 


0.00 ±0.00 


0.00 ±0.00 


Total 


0.96 ±0.32 


0.46 ±0.22 


0.96 ± 0.30 


0.35 ±0.19 


1.48 ±0.48 


1.13 ±0.47 



4 Systematic Uncertainties 

In evaluating systematic uncertainties for the branching fractions, we consider both the uncer- 
tainties that affect the signal efficiency estimate, and uncertainties arising from the background 
estimate. Table |3] lists the systematic uncertainties considered for the signal efficiency. Sources 
of uncertainties that affect the efficiency are: charged-particle tracking efficiency (0.8% per lep- 
ton, 1.4% per charged hadron), charged-particle identification (0.7% per electron pair, 1.9% per 
muon pair, 0.5% per pion), the continuum background and BB suppression selection (1.4%-1.9% 
depending on the mode), and signal simulation statistics (0.1%). The estimated number of BB 
events in our data sample has an uncertainty of 1.1%. An additional systematic uncertainty for 
the efficiency results from the choice of form factor model, which alters the q 2 distribution of the 
signal. We take this uncertainty to be the same as found in [6], which uses the same signal model 
as this analysis. 

Finally, the uncertainties due to the signal efficiency of the tbes an d AE selection requirements 
are determined from the measured mean and width of these distributions in charmonium control 
samples. For B + — > tt + £ + £~ , we use samples of B — > J/ipK + events, in which the mean and width 
are precisely bounded; for B° — ► ir°£ + £~ , we use samples of B — ► J/ipir events, which have limited 
statistics and introduce a total systematic uncertainty of 7%. For the electron modes, we allow for a 
larger or smaller bremsstrahlung tail in the AE distribution, introducing a systematic uncertainty 
of 1-2%. The total systematic uncertainty of the signal efficiency is 5% for B + — > ir + £ + £~ and 9% 
for B° -» 7r°£ + £-. 

Systematic uncertainties for the estimated backgrounds arise from two sources: uncertainties 
of the maximum likelihood fit which determines the combinatorial background, and uncertainties 
of the peaking background estimates. Table El summarizes the sources of background systematics. 
In the maximum likelihood fit to the region defined by 5.2 < tties < 5.2724 GeV/c 2 and \AE\ < 
0.25 GeV, the uncertainty of the extrapolated signal box yield is determined from the induced change 
due to ±lo" variations in the fit region yield, the slope parameter in the tties PDF, and the exponent 
of the AE PDF. We also consider the effect of using different PDF parameterizations on the 
background estimates, and use the computed differences to bound a systematic uncertainty. We fit 
a PDF which is correlated in ?hes and AE via a linear AE dependence in tbes slope parameter. We 
also estimate background from fits for which the AE PDF is a linear or quadratic polynomial. These 
alternative fits result in differences in estimated background generally smaller than the uncertainties 
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Table 3: The sources of systematic uncertainty in signal efficiency (%) considered for ir £ + £ decays. 



Systematic 


7r + e + e 








7r + e/i 


7r u e/i 


Trk eff. 


±3.0 


±1.6 


±3.0 


±1.6 


±3.0 


±1.6 


Electron ID 


±0.7 


±0.7 






±0.4 


±0.4 


Muon ID 






±1.9 


±1.9 


±1.0 


±1.0 


Pion ID 


±0.5 




±0.5 




±0.5 




7T ID 




±3.0 




±3.0 




±3.0 


Fisher and BB likelihood 


±1.4 


±1.4 


±1.7 


±1.9 


±1.4 


±1.4 


MC statistics 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


BB counting 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


Model dep. 


±3.0 


±3.0 


±3.0 


±3.0 


±3.0 


±3.0 


signal toes model 


±0.3 


±5.1 


±0.4 


±4.9 


±0.3 


±5.1 


signal AE model 


±0.6 


±5.1 


±0.5 


±5.4 


±0.5 


±5.2 


signal AE radiative tail 


±1.2 


±1.3 






±1.0 


±1.4 


Total 


±4.9 


±8.8 


±5.1 


±9.0 


±4.9 


±8.9 



of the baseline fit. The uncertainty of the hadronic B decay peaking background is dominated by the 
control sample statistics from which it is derived; the leptonic peaking background uncertainties are 
dominated by the uncertainty of the assumed branching fractions for these processes, particularly 
B+ -^K + e + £~. 

Figures 01 and |1] show projections onto toes and AE of a two-dimensional unbinned maximum- 
likelihood fit to measure the number of background events. The best fit PDFs are consistent with 
the observed data. The fit region yields for each mode are also consistent with the simulated pre- 
dictions from which optimization of selection criteria were derived. We estimate the combinatorial 
background in the signal region from these PDFs and those estimates together with the peaking 
backgrounds result in about one event per decay mode expected in the signal region. 

5 Results 

Figure [5] shows a scatter plot of the fit region with events from the BABAR data that pass all our 
selection criteria. The signal events are expected to be found in the signal region marked by a box 
drawn in the plot. We observe a total of 3 B — > tt£ + £~ candidates and 1 B — > 7re/i candidate in the 
signal region. This is consistent with no significant signal above the expected backgrounds, so we 
calculate upper limits at 90% confidence level via a frequentist method which takes into account 
uncertainties both in the signal sensitivity and in the expected background [19]. The limits are 
presented in Tabled along with all numbers necessary in this method to perform the calculation. 
The upper limits for the electron and electron-muon modes are in the range of 1-2 xlO" 7 ; the muon 
modes are somewhat less sensitive with limits in the range of 2 — 5 x 10 -7 . Assuming the partial 
widths of B — > ir£ + £~ to electrons and muons are equal, the limits for the two decay modes can 
be simply combined to provide a combined limit of 1.06 x 10~ 7 for B + — > ir + £ + £~ and 1.02 x 10~ 7 
for B° — ► ir°£ + £~. Assuming further that there is isospin asymmetry in the partial widths of 
B — > n£ + £~ to charged and neutral pions, the four B — » ir£ + £~ modes can be combined with the 
constraint B(B° —* tt°£ + £~) = \ x t b o/t b + x B(B + — > n + £ + £~) to compute a combined limit, 
expressed in terms of the B + — > ir + £ + £~ branching fraction, of 0.79 x 10 -7 . A similar combined 
limit of 0.98 x 10 -7 is obtained for the lepton-flavor violating mode B — ► tre/j,. 
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Table 4: Results for the B — > ir£ + £~ analysis, including observed signal candidate events, expected 
background, signal yield upper limit at 90% confidence level, signal efficiency, expected branching 
fraction upper limit at 90% confidence level, and the observed branching fraction upper limit at 
90% confidence level. The numbers in parentheses are limits evaluated without the inclusion of 
systematic uncertainties. Combined limits at the bottom are derived from simultaneous limits 
calculated from the individual modes. 







Observed 


Expected 


Events UL 


Signal 


Expected BF UL 


BF UL 




Mode 


Events 


Background 


90% C.L. 


Efficiency 


90% C.L. (KT 7 ) 


90% C.L. (10" 7 ) 


B+ 


— > 7r + e + e _ 


1 


0.96 ± 0.32 


2.97 (2.93) 


7.5 ± 0.4% 


1.64 


1.72 (1.70) 


B° 


-» n°e+e- 





0.46 ±0.22 


1.86 (1.84) 


6.3 ± 0.6% 


1.79 


1.29 (1.27) 


B+ 


— * ir + fi + 


1 


0.96 ±0.30 


2.96 (2.93) 


5.2 ± 0.3% 


2.37 


2.47 (2.45) 


B° 


-> ttV+M" 


1 


0.35 ±0.19 


3.57 (3.55) 


3.4 ± 0.3% 


3.18 


4.56 (4.53) 


B+ 


— > 7r + e^i 


1 


1.48 ±0.48 


2.49 (2.41) 


6.3 ± 0.3% 


2.17 


1.72 (1.66) 


B° 


— ► 7r°e/i 





1.13 ±0.47 


1.28 (1.18) 


3.7 ± 0.3% 


3.52 


1.50 (1.38) 


B+ 


-> n + e + r 












1.06 (1.04) 


B° 


-» *?t+r 












1.02 (1.01) 


B - 


-> n£+r 












0.79 (0.77) 


B - 


-* 7re/i 












0.98 (0.90) 



Table 0] also lists the mean expected upper limits, defined as a sum of upper limits obtained for 
n events observed (n = 0,1,2...), each weighted by the Poisson probability of observing n events 
when expecting background only. 

6 Conclusion 

We present a preliminary result of a search for B — > ir£ + l~ using a sample of (230.1 ± 2.5) x 10 6 
BB pairs produced at the T(4S) resonance. We do not see any excess of events in the signal region, 
and we measure the upper limit at 90% confidence level of the lepton-flavor-averaged branching 
fraction to be 

B(B + -» n + £ + r) = 2 x ^B(B° -» n°£ + r) < 7.9 x 10^ 8 , 

consistent with the Standard Model prediction of 3.3 x 10 -8 . This is the first such search performed 
by the B factory experiments; with anticipated final samples of order 1 ab _1 , and with the small 
backgrounds observed in this analysis, it may be possible in the future to achieve an experimental 
sensitivity comparable to the Standard Model prediction. 
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Figure 3: Projections of mES an d A.E of B —> n£ + £~ events in the full fit region. Superimposed 
is the PDF we use to model combinatorial background. The parameters were obtained from two- 
dimensional unbinned maximum-likelihood fits to events outside the signal region, with mss < 
5.2724 (to the left of the gray vertical line). 
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Figure 4: Projections of mES and AE of B — * irefi events in the full fit region. Superimposed 
is the PDF we use to model combinatorial background. The parameters were obtained from two- 
dimensional unbinned maximum-likelihood fits to events outside the signal region with mES < 
5.2724 (to the left of the gray vertical line). 
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Figure 5: Distribution in two dimensions (mES and A.E) of events in the background fit region 
passing the selection. The small box defines the signal region from which signal candidates are 
selected. 
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